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Claims 



10 



1. A method for producing three-dimensional information of an object (4) in 
medical X-ray Imaging, characterized in that 

- the object is modelled mathematically independentiy of X-ray imaging, 

- the object is X-radiated from at least two different directions and the said 
X-radiatlon is detected to form projection data of the object (4), 

- and said projection data and raid mathematical modelling of the object are 
utilized In Bayesian inversion based on Bayes' formula 



Ppr ix)P(jn I 



to produce three-dimensional fnfonnatlon of the object;, the prior distribution 
ppr(x) representing mathematical modelling of the object, x representing the 
object image vector, which comprises values of the X-ray a^nuation 
coefficient Inside the object, m representing projection data, the likelihood 
15 distribution p(m|x) representing the X-radiatlon attenuation model between 
the object image vector x and projection data m, p(m) being a normallzatioh 
constant and the posteriori distribution p(xlm) representing the three- 
dimensional Information of the object (4). 

20 2. A method aax>rxllng to dalm 1, characterized in that the three- 
dimensional information of the object (4) is one or more two-dimensional 
Images representing X-ray attenuation coefficient along slices through the 
object. 

25 3. A method according to claim 1, characterized in that the tiiree- 
dimensional information of the object (4) is a three-dimensional voxel 
representation of the X-ray attenuation in the object 



30 



4. A method according to claim 1, characterized in that the measurement 
model is m = Ax+e, where matrix A contains tt>e lengths of the path of the 
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X-ray inside each voxel and the noise e Is independent of object image 
vector X leading to the likelihood distribution 



5 5. A method according to dalm 1, characterized in that the said 

mathematical modelling comprises that X-radiation attenuates when passing 
the object (4), which means that every image voxel Is nonnegative. 

6. A method according to claim 1, characterized In that mathematical 
10 modelling Is expressed by the formula: 

Ppr(a?) = exp(-a ^ Uf^(x)) 

N 

Where the sum is taken over a collection of 3D neighbourhoods N and the 
15 value ^« depends only on the values of voxels belonging to the 

neighborhood N, and a is a positive regularization parameter used to tune 
the width of the prior distribution. 

7. A method according tD claim 1, characterized In that the 3D 

20 tomographic problem Is divided Into a stack of 2D tomographic problems and 
on every such 2D problem, the mathematical modelling is expressed by the 
formula: 

Ppr(«) = exp(-a ^ C/jv (a?)) 

N 

25 where the sum Is taken over a collection of 2D neighbourhoods N and the 
value depends only on the values of pixels belonging to the 
neighbortiood and a Is a positive reguIarizaOon parameter used to tune 
the width of the prior distribution, and the 2D tomographic problems are 
related to each other by the formula 

30 

pr3D(xa)) = exp(- y ZS I xO)[Kq] - xQ-l)[Kq] \), 
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where the sums are taken over all pixels (k=l,...,l^ q=l,...,Q ) and y > 0 Is 
another regularization parameter. 

8. A method according to claim 7, characterized in that the neighborhoods 
5 consist of two adjacent pixels and U calculates a power of the absolute value 

of the difference, leading to the fonmula 

10 +i:i:i^<'*>lfe,9lT:.a^<'*nfc,9 + l]r)) 

k=l q=l 

where s is a positive real number. 

9. A method according to claim 8, characterized in that s=l corresponding 
15 to tx)tal variation (TV) prior desalbing objects (4) consisting of different 

regions with well-defined boundaries. 

10. A method according to daim 1, diaracterized In that mathematical 
modelling is qualitative structural information of the target where the 

20 structural information is encoded In prior distributions that are concentrated 
around object image vectors x that correspond to the physiological structures 
of the object (4). 

11. A method according to daim 1, characterized in that mathematical 
25 modelling consists of a list or probability distribution of possible attenuation 

coeffldent values in the object (4). 

12. A method according to daim 1, characterized in that the X-ray 
imaging geometry, such as X-ray source position, has unknown error 

30 modelled in the distribution p(mix). 
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13. A method according to daim 1, characterized in that the X-radiation 
measurement noise is Poisson distributed. 

14. A method according to daim 1, characterized in that the medical X-ray 
5 imaging is dental radiography. 

15. A method according to daim 1, characterized in that the medical X-ray 

,10 .16. A method according to daim 1, characterized in that the medical X-ray 
imaging is mammography. 

17. A method according ta daim 1, characterized in that three- 
dimensional Information of the object (4) Is produced oh the basis of the 
15 maximum a posteriori estimator (MAP) which is calculated by the equation: 



m representing projection data and x representing the object Image vector 
and where the maximum on the right hand side of the equation Is taken over 
20 all X. 

18. A method according to daim 1, characterized in that three- 
dimensional Infomrration of the object (4) is produced on the basis of the 
conditional mean estimator (O^), which Is calculated by the equation: 



p(iPMAP I m) — maxp(x \ m), 



25 




where m represents projection data and x represents the object image 
vector. 



30 
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19. A medical X-ray device (5) arrangement for producing three-dimensional 
Informaf on of an object (4) in a medical X-ray imaging, characterized in 
that the medical X-ray device (5) arrangement comprises: 
- means (15) for modelling the object (4) mathematically Independently of 



- an X-ray source (2) for X-radlating the object from at least two different 
directions 

- a detector (6) for detecting the X-radlatlon to form projection data of the 
ob|ect(4) 

10 - and means (15) for utilizing said projection data and said mathematical 
modelling of tiie object In Bayesian inversion based on Bayes' formula 



15 to produ(% three-dimensional Infonfnation of the object, the prior distribution 
Pp,(x) representing mathematical modelling of the object, x representing tiie 
objed: Image vector, which comprises values of the X-ray attenuation 
coeflRdent Inside the object, m representing projection data, tiie likelihood 
distribution p(m|x) representing the X-radiation attenuation model between 

20 the object image vector x and projection data m, p(m) being a normalization 
constant and the posteriori distribution p(xim) representing the three- 
dimensional infonnation of the object (4). 

20. A medical x-ray device (5) arrangement according to dalm 19, 
25 characterized in that the three-dimensional Information of the object (4) is 
one or more two-dimensional Images representing X-ray attenuation 
coeffident along slices through the object. 



5 



X-ray Imaging 



Ppr(x)p(m\x) 



30 



21. A medical x-ray device (5) an^ngement acconding to daim 19, 
characterized in that the three-dimensional Information of the object (4) is 
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a three-dimensional voxel representation of the X-ray attenuation in the 
object 



22. A medical X-ray device (5) arrangement according to claim 19, 

5 characterized in that the medical X-ray device arrangement comprises 
means (15) for modelling the mea^rement as 
m = Ax+e, 

where matrix A contains the lengtlis.of the path of the X-ray Inside each 
voxel and the noise e is indiependent of object Image vector x leading tn the 
10 likelihood distribution 

m 

23. A medical X-ray device (5) arrangement according to claim 19 

15 characterized in tiiat the medical X-ray device arrangement comprises 
means (15) for modelling the object (4) mathematically so that X-radlat*on 
attenuates when passing the object (4), which means that every Image voxel 
is nonnegative. 

20 24. A medical X-ray device (5) an^ngement accoreling to claim 19, 

characterized in that the medical X-ray device arranganent comprises 
means (15) for modelling the object (4) mathematically by the formula: 

Ppr(a?) = exp(-a Un{x)) 
25 JV 

where the sum is taken over a collecBon of 3D neighbourhoods N and the 
value ^"^^^ depends only on the values of voxels belonging to the 
neighborhood N, and a Is a positive regularization parameter used to tune 
the width of the prior distribubon. 

30 

25. A medical x-ray device (5) an^ngement according to claim 19, 
characterized in that the 3D tomographic problem is divided into a stack 
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Of 2D tomographic problems and on every such 2D problem, and the medical 
X-ray device arrangement comprises means (15) for modelling the object (4) 
matJiemab'cally by the formula: 

5 Ppr(») = exp(-Qj ^ Un {x)) 

N 

Where the sum is taken over a collection of 2D neighbourhoods N and the 

value W depends only on the values of pixels belonging to the 
neighborhood and a is a positive regularizaOon parameter used to tune 
10 the width of the prior distribution; and the 2D tomographic problems are 
related to each other by the formula 

pr3D(xa)) = exp(- y ZS I xa)[k,q] - xa-l)[k,q] I)/ 

15 where the sums are taken over all pixels (k=l,...,K, q=l,...,Q ) and y > 0 is 
anotiier regulatlzadon param^sr. 

26. A medical X-ray device (5) arrangement according to daim 25, 
characterized in that the neighboitiood systems oon^t of two 

20 neighboring pixels xj, )* or voxels xj, xk and calculates a power of the 

+EE ^"'ftsi -«<»t*.9+iin) 

25 

absolute value of the difference, leading to the formula where s is a positive 
real number and a Is a regularization parameter used to tune the width of 
the prior distribution. 

30 27. A medical X-ray device (5) arrangement according to claim 26, 

ciharacterized in that the medical X-ray device arrangement comprises 
means (15) for modelling the object (4) mathematically by setting s=l 
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corresponding to total variation (TV) prior describing objects consisting of 
different regions with well-defined boundaries. 

28. A medical X-ray device (5) arrangement acconding to daim 19, 

5 characterized in that the medical X-ray device arrangement comprises 
means (15) for modelling the object (4) mathematically by assuming that 
mathematical modelling is qualitative structural information of the target 
where tlie structgral information is encoded In prior distributions that are 
concehtrated around image vectors x that correspond to the physiological 
10 structures of tiie target 

29. A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray device arrangement comprises 
means (15) for modelling the object (4) mathematically by assuming that 

15 mathematical modelling consists of a list of possible attenuation coefficient 
values in the object. 

30. A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray device arrangement comprises 

20 means (15) for modelling the object (4) mathematically by assuming tiiat the 
X-ray imaging geometry, such as X-ray source position, has unknown error 
modelled in the distribution p(m|x).. . 

31. A medical X-ray device (5) an^ngement according to daim 19, 

25 characterized in that the medical X-ray device arrangement comprises 
means (15) for modelling the object (4) mathematically by assuming that X- 
radiation measurement noise is Poisson distributed. 

32. A medical X-ray device (5) an^ngement according to daim 19, 

30 characterized In that the medical X-ray imaging is dental radiography. 
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33. A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray imaging is surgical C-arm imaging. 

34. A medical X-ray device (5) arrangement according to claim 19, 

5 characterized in that the medical X-ray imaging is mammography. 

35. A medical X-ray device (5) arrangement according to claim Id, 

: characterized in that the medical X-ray device arrangement comprises 
means (15) for producing three-dimensional infonmatlori of the object (4) on 
10 the basis of the maximum a posteriori estimator (MAP), which is calculated 
by the equation: 



15 m representing projection data and x representing the object image vector 
and where the maximum on the right hand side of the equation is taken over 



36. A medical X-ray device (5) arrangement according to claim 19, 
20 characterized in that ttie medical X-ray device arrangement comprises 
means (15) for producing three-dimensional information of the object (4) on 
the basis of the conditional mean estimator (CM), which is calculated by the 
equation 



p(iPMAP I 'nri) = maxp(x \ m), 



all X. 



25 



XcM = J xp(x\m)dx 

where m represents projection data and x represents IJie object image 
vector. 




